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The correct quantum description for curvature squared term in the action can be obtained by 
casting the action in the canonical form with the introduction of a variable which is the negative 
of the first derivative of the field variable appearing in the action, only after removing the total 
derivative terms from the action. We present the Wheeler-deWitt equation and obtain the expression 
for the probability density and current density from the equation of continuity, further in the weak 
energy limit we obtain the classical Einstein's equation. Finally we present a solution of the wave 
equation. 
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O ■ I- INTRODUCTION 

^ . 

In analogy with the particle quantum dynamics and the quantum field theory,Hartle-Hawking EJ suggested that 
the ground state wave function of the universe is given by 



Mhij}=N 5gex P {-I E [g]), (1) 



O 
O 

with the proposal that the sum should be taken over all compact four geometries. This implies that the universe 
does not have any boundary in the cuclidean regime. The above functional integral over all compact four geometries 
bounded by a given three geometry is interpreted as amplitude of the given three geometry to arise from a zero three 
geometry, where ipo[hij] is finite. Thus the ground state is the amplitude of the universe to arise from nothing. 

The euclidean form of the Einstein-Hilbert action is not positive definite. Therefore the functional integral for the 
j^ , wave function of the universe as proposed by Hartle-Hawking runs into serious problems, since the wave function 
&Q' diverges badly.There are some prescriptions to make the above path integral convergent. In one of the prescriptions 
@i Hi; (§' ft was suggested that the proper choice of the contour of integration in the functional integral might 
yield positive definite action, while in another prescription fq] a suitable choice of conformal factor may give rise to 
positive definite action. Unfortunately, a completely satisfactory result has not yet been obtained through the above 
two prescriptions. 

Further to explain the inflationary scenario without invoking phase transition in the very early universe, Starobinsky 
P| considered a field equation containing only the geometric terms, but this field equation could not be obtained from 
the action principle, as the terms in the field equations are obtained from the perturbative quantum field theory. 
Starobinsky and Schmidt pj] later have shown that the inflationary phase is an attractor in the neighbourhood of the 
solution of the fourth order gravity theory. 

However, to get a convergent integral for the wave function Horowitz S proposed an action in the form 



S = - J d 4 XV=g~{AC% kl +B(R- 4A) 2 ] , (2) 



where, Cijki is the Weyl tensor, R is the Ricci scalar, A is the cosmological constant and A, B are coupling constants. 
In the weak energy limit this action (2) reduces to the Einstein-Hilbert action for gravity. To obtain a workable and 
simplified form of the field equations, we consider a homogeneous and isotropic minisuperspace model. Then Ciju 
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vanishes trivially and if one chooses A = 0, for the sake of simplicity, the action (2) retains only the curvature squared 
term. 

The classical held equations for R 2 gravity can be obtained by the standered variational principle. In the variational 
principle total derivative terms in action are extracted and yields a surface integral and it is assumed that the surface 
integral vanishes at boundary of variation, or, the action is chosen in such a way that those surface integral terms 
have no contribution in the action. However for canonical quantisation this principle is not of much help. For this 
purpose one has to express the action in the canonical form. Boulware et al B] has proposed that this can be achived 
by introducing a new variable, which is the hrst derivative of the action with respect to the highest derivative of the 
held variable that appears in the action. 

Hawking and Luttrell |L(| utlised the technique to identify the new variable and the conformal factor and showed that 
the Einstcin-Hilbert action along with the curvature squared term reduces to the Einstein's gravity being coupled 
to a massive scalar held. Horowitz on the other hand retained only the curvature squared term and showed 
that the corresponding Wheeler-deWitt equation looks very similar to the Schrodinger equation. The semi-classical 
approximation to such quantum equation has also been presented by him. Pollock |ll|]has also introduced the same 
technique to the induced gravity theory and obtained the same type of result as that obtained by Horowitz Q in the 
sense that the corresponding Wheeler-deWitt equation looks similar to the Schrodinger equation. 

In the present paper we shall deal with the curvature squared term with a view to study at which stage the 
new variable should be introduced to obtain the correct Wheeler-deWitt equation. Our proposal is that one should 
introduce the new variable only after removing all the total derivative terms from the action and the new variable is 
the first derivative of the action (free from total derivative terms) with respect to the highest derivative of the field 
variable. Then one should express the action in the canonical form with respect to the new variable. At this stage 
one should also remove the total derivative terms (if any) from the action and the variation of the action yields the 
classical field equations. 

It is important to note that the classical field equations remain the same both in our proposal and in the Boulware's 
H proposal. However the canonical quantization in our proposal yields inequivalent Wheeler-deWitt equation from 
that obtained by the technique followed by Horowitz || . 

To elucidate our proposal we shall consider a toy model of the vacuum Einstcin-Hilbert action in section 2., and show 
that the introduction of a new variable following Horowitz and Pollock leads to a wrong Wheeler-deWitt equation. 
We therefore conclude that such a variable should be introduced only after removing all the total derivative terms 
from the action. 

In section 3, we shall consider curvature squared action and with the proper choice of the variable ,we show that 
the corresponding Wheeler-deWitt equation is different from that obtained by Horowitz ||, though it still looks 
similar to the Schrodinger equation. In order to compare these results we have chosen FRW model in conformal time 
coordinate for closed (k = 1) universe. The Wheeler-deWitt equation thus obtained contains an effective potential V. 
The presence of the potential is in sharp contrast from that obtained by Horowitz . The potential V is a function of 
the expansion parameter (i.e. the Hubble parameter) of the classical field equations and it is diverging both at the 
zero's and at the very large values of the expansion parameter, further the effective potential is also asymmetric with 
respect to the expansion and collapse of the universe. The physical interpretation of the potential is not yet clear to 
us, however its extremum yields interesting results. One of the extremum of the potential corresponds to the classical 
vacuum Einstein equation , while the other can be viewed upon as a classical equation which admits inflation . In fact, 
the classical field equation obtained from the extremum of the potential is not surprising, but it should be a desirable 
feature of higher derivative gravity theory at the weak energy limit (i.e. at the long wavelength limit), where the 
contribution of the kinetic energy is sufficiently small with comparison to the potential energy in the Hamiltonian; 
justifying the idea of extremum of the potential. 

In section 4, we shall present the Hamilton-Jacobi equation corresponding to the Schrodinger equation and then we 
can identify the time parameter that appears in the Hamilton-Jacobi equation. In this section we shall also present 
a particular solution of the Wheeler-deWitt equation for n = operator ordering, which is well behaved for all values 
of the expansion parameters and at the vanishing three space volume. The probability density ip*ip is a real quantity 
and is proportional to the scale factor of the universe. 

In section 5, we have treated the same problem in proper time coordinate for arbitrary three space curvature 
(k = o, ±1) . This shows that the form of the Wheeler-deWitt equation remains unchanged . We then turn our 
attention (see section 6) to compare our formalism to that developed by Schmidt [p.2| earlier. For this purpose we 
have chosen k = and have shown that the formalism of Schmidt |12| following Ostrogradski [reference cited therein] 
is different from that of ours, though the new variables in both the formalisms turn out to be the same and hence 
the Hamiltonian. This approach of Schmidt has recently been considered by Fabris and Reuter J13| . However in our 
method we have introduced the variables in two stages systematically first, to obtain the field equations and then to 
find the right Hamiltonian that should be quantized. The method is simple and straight forward. While the other 
formalism introduced the variables at one stage only, to find the Hamiltonian. This formalism would yield the same 



Hamiltonian as that of ours only if one starts from an action free from surface term. 
In section 7, we present a brief conclusion. 

II. INTRODUCING THE NEW VARIABLE IN THE VACUUM EINSTEIN-HILBERT ACTION (A TOY 

MODEL) 

The vacuum Einstein-Hilbcrt action is 

1 



S = -^J ^X^—gR. (3) 

We consider the closed Robertson- Walker metric in the form 

ds 2 = exp2a(r))[dr] 2 - d X 2 - sin 2 X (d9 2 + sin 2 dd(f> 2 )}, (4) 

da 
dq ' 



where R = — 6exp[— 2a] (1 + a 2 + a) and a = ^l. 



The action (3) can be expressed as 

S= exp[2a](l + a + a)dri = / exp[2a](l — a )drj + S mi (5) 

where S m is the surface term given by S m — exp[2a]a.The classical field equation (assuming S m = at the 
boundary) is 

a + a 2 + l = Q. (6) 

The constraint equation is H = Oi.e., 

1 + a 2 = (7) 

and the Wheeler-deWitt equation is 

d 2 ip ndip 4 IA , 

d^ + alk;-]? 6 ^^ ' (8) 

where, n is the operator ordering index. Now instead of removing total derivative term from the action (5) if one 
introduces a new variable Q as 

then we have Q = —exp[2a\. Hence the action (5) can be written in the canonical form as 

S=- [ [Q(l + a 2 ) + Qa]d v = f [Qa - Q(l + a 2 )]d v - S m . (10) 

one can observe that the surface terms both in (5) and (10) are the same. Apart from the surface term in (10) we 
have 

dL 
p a = ^T=Q-2Qa (11) 

and 

oQ 
Then the classical field equation with respect to the variation of a is 

Q-2Qa-2Qa^0. (13) 



One can observe that the above field equation is identically satisfied by equation Q = —exp(2a). The other field 
equation is 

a + a 2 + 1 = 0. (14) 

Hence, the classical field equations do not change even with the introduction of the new variable. We now express 
the Hamiltonian constraint in the phase space variables, so 

H = Pa p Q + Qpq + Q = 0. (15) 

Now canonical quantization should be performed with basic variables viz, a and a (Boulware et al). Let a = .t,so 

~jS _ _9L _ PL 
a da dx 

H = xp a - x 2 p x - p x = 0. (16) 



we replace pq by x and since Q = — §§ = — §§ = ~ fj = ~Vx , so Q should be replaced by p x . Hence 



Hence the Wheeler-deWitt equation in this case turns out to be 

dip 1 + x 2 dtp 



(17) 



da x dx ' 

whose solution is 

ip = ip exp(ca){l + x 2 ) c/2 , (18) 

where c is a constant. Clearly this does not satisfy equation (8). Hence the equation (17) is not the correct 
Whcclcr-dcWitt equation. Therefore we conclude that in order to obtain the correct quantum dynamics one should 
first remove all the total derivative terms from the action and then introduce the new variable if at all required. 

III. CURVATURE SQUARED ACTION AND THE CORRECT WHEELER-DEWITT EQUATION 

We assume the action as 



S = -— {l + a' + aydr], (19) 

where the background geometry is given by metric (4). According to our proposal we remove the total derivative 
terms from the action (19).Hence 

S =- 1 jf [(1 + « 2 ) 2 + u 2 ]d v + S m7 (20) 

where m = 12n 2 (3 and S m = — Tf (3 + a 2 ) .Now we introduce the new variable to cast the action (20) in the canonical 
form 

^ = -g=ma/2. (21) 

Hence 

S = m/AJ [aQ-{l + a 2 ) 2 + Q 2 /A]df 1 + S mi (22) 

where S mi is the surface term given by 

S mi = -— d(3 + a 2 + 3d). (23) 

6 

Assuming S mi vanishes at the boundary, the classical field equations are 

Q - 12d 2 d - 4d = 0, (24) 

and 



Q = 2d. (25) 

The definition of Q in (21) is thus recovered in (25) and the equation (24) is the correct classical field equation as 
can be observed by replacing Q by 2d . Now the Hamiltonian constraint equation is 

4 64 a 8 9 to ^ 9 to . _ 

H=- PaPQ + — p 4 + £ Q 2 + -r, 26 

TO TO d V TO ^ 16 4 

where p Q and pg are the canonical momenta corresponding to a and Q respectively. Since the canonical quantization 
should be performed with basic variables viz., a and d, therefore choosing a — x, we replace pq by ^f- and Q by 



•^2a .Hence 



F = xp Q -^ + ^(l + x 2 ) 2 . (27) 



to 4 



The corresponding Wheeler-deWitt equation is 



. % dip h d 2 ijj nh dip tn . 2 \2 (r, Q \ 

da mx dx 2 mx 2 dx Ax 

where n is the operator ordering index. Though the Wheeler-deWitt equation (28) looks similar to the Schrodinger 
equation it differs from that obtained by Horowitz |g). For n = 0, the first derivative term in the right hand side of 
(28) disappears, while the Wheeler-deWitt equation obtained by Horowitz always retains the first derivative term on 
the right hand side. We have also obtained an effective potential V(x) in the process, which has not been obtained by 
Horowitz. The equation (28) can be written as 

ih^ = Hoi;, (29) 

where, 



h 2 , d 2 n d 



in 



H o = — (7T^ + -7r) + ^( 1 + x ) ■ (30) 

mx ox x ox 4x 

It is noted that Hq operator is hermitian and as a consequence we obtain the equation of continuity 

!^+V-J = 0, (31) 

where p and J are the probability density and the current density respectively for n = — l.The p and J are given 
by p = ip*tp and J = (J x , 0, 0); where 

Jx = — W>V*-tWS)- (32) 

mx 

For other operator ordering index n we can also obtain the continuity equation, but with respect to a new variable 
y which is functionally related to x only. It is noted that no such continuity equation can be obtained from the 
Wheeler-deWitt equation presented by Horowitz || . In analogy with the quantum mechanics it is apparent from (29) 
and (31) that the variable a can be identified as the time parameter in quantum comosmology, and the variable x acts 
as a spatial coordinate variable in the quantum cosmology. The coordinate variable x = d, which is just the expansion 
parameter and the Hamiltonian operator Hq diverges at the bounce of the universe i.e., d = 0. Further in the weak 
energy limit, the contribution of the kinetic energy is sufficiently small with respect to the potential energy in the 
Hamiltonian Hq. At that regime the Hamiltonian is almost dominated by the potential energy. Now the extremum 
of the potential V(x) yields 

1 + d 2 = (33) 

or, 

d 2 - 1/3 = 0. (34) 

The equation (33) is the classical vacuum Einstein equation. It is to be noted that the equation (33) which we obtain 
here is not surprising, rather it should be a desirable feature of higher derivative gravity theory at the weak energy 
limit. Equation (34) on the other hand is just the equation admitting inflation, supporting the idea of Starobinsky to 
explain the inflationary scenario without invoking phase transition. 



IV. SOLUTION OF THE WHEELER-DEWITT EQUATION 

The Hamilton Jacobi equation can be obtained from (29) and is given by 

dA 1 ,dA. 2 m, 2 . 2 . ,„_ s 

-7T + — -~-) 2 -— l + a: 2 ) 2 = 0, 35 

da rax ox Ax 

where A is the Hamilton Jacobi function. This is the exact form of the classical Hamiltonian equation given by 

dA 

^+^ = 0, (36) 

2 

where H is the Hamiltonian and is given by H = -^*- — j^(l + x 2 ) 2 . The parameter a in (36) or (35) is simply 
the time parameter as we have already claimed from equations (28) and (29) in the previous section. 
A particular solution of (28) for n — is given by 

%Tfl X 

ip = ip exp[a + —(x + — )]. (37) 

The wave function ip is well behaved for all values of x (i.e., for any expansion parameter) and at the vanishing 
three space volume. The probability density ip*ip is just \i(>o\ 2 exp(2a) , and it is a real quantity and propor tional to 
the scale factor of the universe. However for normalization of the probability density one needs boundary condition 
of the universe. Recntly, Bachmann and Schmidt |l4[ have solved the Wheeler-deWitt equation for a minimally 
coupled scalar held in Bianchi type-I model and shown that the notion of a cosmological quantum boundary is well 
defined. The solution (37) is not a semiclassical one, as the terms within the exponential of (37) does not satisfy 
the Hamilton-Jacobi equation. Therefore it is not clear wheather such boundary can be called upon in the present 
context. 

V. CURVATURE SQUARED ACTION IN PROPER TIME COORDINATE 

In this section our first attempt is to show that the Wheeler-deWitt equation remains unchanged even in the proper 
time coordinate and then we compare our formalism with that of Schmidt [ |12[ formulated earlier. For this purpose 
we consider the FRW metric in proper time coordinate given by 

dv 

ds 2 = dt 2 - a 2 \- — - + r 2 (d8 2 + sin 2 6d(t) 2 )}. (38) 

1 — kr z 

Removing the total derivative terms the action can now be cast in the following form 

S = --[ [aa 2 + id2 + k)2 ]dt + S 1 , (39) 

A J a 

Where dot represents derivative with respect to the proper time and Si = — ^(fca + d 3 /3). Introducing the new 
variable 

Q = --t? = 2 ^ (40) 

m oa 

and expressing the action in the canonical form we have 

5 = ?/ [aQ-^±^ + ^ ]d t + S 2 , (41) 

where S2 — — ^d(3fc + a + Sad). The canonical momenta are 

V*= d -^r = ^[Q- Ma 2 + k)/a] (42) 

oa 4 



and 



9L 

PQ = ^X = "W 4 - (43) 

oQ 

It is now possible to find the classical field equations, one of which turns out to be the definition of Q and the other 
is found to resemble with one already obtained in section 3. The Hamiltonian turns out to be 



4 . l6p 4 Q 2kp 2 m 2 k 2 m 2 Q 2 

— {PaPQ + —j- + + -T-„ -^7— 

m m z a a Iba 64a 



ff =>^ + ^ + — + ! 7^-^)- ( 44 ) 



Now, for canonical quantization we invoke the basic variables, i.e., choose a = x, replace Q by —Ap y /m and pq by 
my/4. Thus the Hamiltonian turns out to be 

H = yp y -^ + ^(y 2 + k) 2 , (45) 

ma 4a 

whose quantum version viz., the Wheeler-deWitt equation is 

oa my oy z y ay Ay 
where n is the operator ordering index. Now choosing a — e q the Wheeler-deWitt equation can be written as 

..dib h 2 ,d 2 ib ndth. m, 2 ,si , ,„, 

oq my oy z y ay Ay 

It is interesting to note that this form of Wheeler-deWitt equation is the same as that obtained in equation (28) of 
section 3, in conformal time. 

VI. COMPARISON WITH SCHMIDT'S FORMALISM 

Let us turn our attention to compare our formalism with that proposed by Schmidt |12J earlier. For this purpose 
we choose k — and a — e 9 .Thus the action (39) turns out to be 

S = -j f fe^dt. (48) 

Following the method already proposed, our new variable is 

Q = 2qe 3q (49) 

and the canonical action is 

711 f r,.A , <9 2 _-3 31 



S =J [qQ + ^e- iq }dt + S 3 , (50) 

vwhere S3 = —^-Qq- The canonical momenta are p q = mQ/A and pq = mq/A. Hence the Hamiltonian is 

H = - PqPQ - ^Q 2 e~ 3 ". (51) 

m 16 



Introducing basic variable x = q and replacing Q by —Ap x /m and pq by mx/A, we find 

m 



H = x Pq - ^e~ 3q . (52) 



Now replacing qbyQi,q = x by Q 2 , p q = mQ/A = ff t (2qe 3q ) by -mP l /A and finally p x = -mQ/A = -f(2qe 3q ) 
by —mPijA we can express equation (52) as 

H = -^{ Pl Q 2 + \p 2 e- 3 ^). (53) 



This Hamiltonian other than the factor — m/4 was obtained by Schmidt |12] following Ostrogradski's prescription. 
Though the new variables introduced in both the cases to construct the Hamiltonian are found to match, yet the 
formalisms differ. In our method after eliminating surface terms from the action we have introduced variables at two 
stages, first to cast the action in the canonical form for obtaining the classical field equations and then to construct the 
Hamiltonian suitable for quantization, basic variables have been called upon. On the other hand the basic variables 
in the Schmidt's proposal have been introduced from the very begining and there is no prescription suggested to 
eliminate the surface term. Here lies the difference between the two proposals. The Hamiltonian (53) turns out to be 
the same since Schmidt has considered the action (48) which is already free from the surface term. It is not difficult 
to show that Schmidt's formalism would yield a different Hamiltonian if action (19) would have been considered. 

VII. CONCLUSION 

We conclude that whenever a Lagrangian contains terms higher than the first derivative, the correct quantum 
dynamics would be obtained if one casts the action in the canonical form by introducing a new variable only after 
removing all the total derivative terms from the action. We have shown that in our proposal there is no change in 
the classical field equations, but the quantum dynamics changes sharply as it is evident from (29) in comparison with 
the same equation obtained by Horowitz. Therefor, the quantum dynamics is greatly influenced by the presence of 
the surface term in the action. Further the equation of continuity (31) identifies the nature of space and time like 
variables in the minisuperspace and therefore establishes the idea of probality density and current density in quantum 
cosmology. In the weak energy limit we are able to recover the classical Einstein's field equations from the quantum 
dynamics of higher derivative gravity theory. The solution (37) corresponds to a real probability density and it is not 
a semiclassical solution, however for normalization one needs boundary conditions of the universe. 
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